Experimental verification of PbBi 2 Te 4 as a 3D topological insulator 
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The first experimental evidence is presented of the topological insulator state in PbBi2Te4. A 
single surface Dirac cone is observed by angle-resolved photoemission spectroscopy (ARPES) with 
synchrotron radiation. Topological invariants Z2 are calculated from the ab initio band structure to 
be 1; (111). The observed two-dimensional iso-energy contours in the bulk energy gap are found to 
be the largest among the known three-dimensional topological insulators. This opens a pathway to 
achieving a sufficiently large spin current density in future spintronic devices. 

PACS numbers: 73.20.-r, 79.60.-i 



Topological insulators (TIs) have recently emerged as 
a new state of quantum matter, which are distinguished 
from conventional insulators by having a surface state 
with a massless Dirac dispersion in the bulk energy gap. 
The spin orientation of the topological surface state is 
locked with its crystal momentum, resulting in a helical 
spin texture. This new state can be classified by so-called 
Z 2 topological invariants [j], Q ■ Generally, the massless 
Dirac cone can be created at an interface of two materi- 
als: a topological and an ordinary insulator. The unique 
properties of topological surface electrons provide a fer- 
tile ground to realize new electronic phenomena, such as a 
magnetic monopolc arising from the topological magneto- 
electric effect and Majorana fermions at the interface 
with a superconductor 0, 0). Owing to time-reversal 
symmetry, topological surface states are protected from 
backscattering in the presence of a weak perturbation, 
which is required for the realization of dissipationless spin 
transport without external magnetic fields in novel quan- 
tum devices @] . 

A number of materials that hold spin-polarized sur- 
face Dirac cones have been intensively studied, such 
as Bi^Sb, Bi 2 Te 3 iEl, Bi 2 Se 3 EMS, and 

thallium-based compounds [15l - l2l| . Among them, Bi 2 Sc3 
has been regarded as one of the most promising candi- 
dates for potential applications in ultra-low power con- 
sumption quantum devices that can work stably at room 
temperature due to a sufficiently large bulk energy gap. 
Although significant efforts have been made towards 



spintronic applications, the surface contribution to con- 
duction was hardly observed even at low bulk carrier den- 
This motivates the search of new topologi- 



22 24 



sity 

cal insulators with higher density of spin-polarized Dirac 
fermions. 

Recently, some of the Pb-based ternary chalcogenides 
have been proposed as 3D topological insulators [251427 1. 
The present study is focused on PbBi2Te4, which has 
the simplest crystal structure among them. The crystal 
is composed of seven-layer blocks with the atomic layer 
sequence Tc-Bi-Tc-Pb-Te-Bi-Te [sec Fig. SI (a) of sup- 
plementary material (SM)] [2^ |. The theoretical analy- 
sis of Ref. [22| indicates that this compound is a three- 
dimensional (3D) topological insulator. In this Letter, we 
report the first experimental evidence of the topological 
surface state in PbBi 2 Te4 by angle-resolved photoemis- 
sion spectroscopy (ARPES). We establish that the size of 
the two-dimensional iso-energy contour, which is propor- 
tional to the surface carrier density, is the largest among 
the known 3D topological insulators. Thus, PbBi 2 Tc4 
can be counted as one of the most promising candidates 
for realizing a large spin current density in future spin- 
tronic devices. 

A single crystalline sample of PbBi 2 Te4 was grown by 
the standard procedure using Bridgman method [see SM 
for more details] . For ARPES measurement the samples 
were in situ cleaved along the basal plane. Photoemis- 
sion experiment was performed with synchrotron radi- 
ation at the linear undulator beamline (BL1) and the 
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FIG. 1: (Color online) (a) ARPES energy dispersion curves along the TK line acquired at hv = 7.5 eV. Constant energy contour 
maps (i) and their second derivatives (ii) obtained from the ARPES measurement at hv = 7.5 eV (b) and 10 eV (c) in the kii 
range —0.3 A -1 < k x ,k v < +0.3 A -1 from Eb = 470 meV (Dirac point) to (Fermi level). The three-fold symmetrization 
procedure is applied only for the second derivative images, (c) Three dimensional map for hv = 10 eV. 



helical undulator beamline (BL9A) of Hiroshima Syn- 
chrotron Radiation Center (HiSOR). The ARPES spec- 
tra were acquired with a hemispherical photoelectron an- 
alyzer (VG-SCIENTA R4000) at 17 K. The overall energy 
and angular resolutions were set to 10-20 meV and 0.3°, 
respectively. 

Figure 1(a) shows the ARPES energy dispersion curve 
along the f K line of the surface Brillouin zone [see SM, 
Fig. Sl(b)] measured at hv = 7.5 eV. The parabolic 
band with the energy minimum at the binding energy 
of Eb = 200 meV is seen to exhibit a rather strong pho- 
toemission intensity near the T point. The position of 
the energy minimum docs not change with the photon 
energy [see SM, Figs. S2(a) and S2(b)], which confirms 
the two-dimensional nature of this electronic state. More 
importantly a linearly dispersing feature, i.e. the Dirac 
cone with the crossing point at Eb = 470 meV is ob- 
served. 

Figures 1(b) and 1(c) summarize the constant energy 
contours (i) and their second derivatives (ii) in the ky 
range -0.3 A^ 1 < k x ,k y < +0.3 A" 1 from E B = 
470 meV (Dirac point) to (Fermi level) at hv = 7.5 
and 10 eV, respectively. With hv = 10 eV we find at 
the Dirac point energy six ellipses oriented along f M in 
addition to the point-like feature at the f point for the 
Dirac state. On the other hand, the elliptical contours 
are much weaker at hv = 7.5 eV, signifying a strong ma- 
trix elements effect. In going away from the Dirac point, 



the sole hexagonally shaped contour is observed down to 
Eb = 290 meV and at smaller binding energies another 
state becomes enclosed inside the Dirac cone (starting 
with Eb = 260 meV). It is interpreted as the onset of the 
bulk conduction band (BCB) because the shape of the 
inner contours considerably depends on the photon en- 
ergy, see Figs. 1(b) and 1(c). Below Eb = 160 meV the 
Dirac cone further deforms and practically merges into 
the bulk conduction band. In the map measured with 
hv = 10 eV the hexagonal Fermi surface of the Dirac 
cone encloses two large and one small triangular shaped 
surfaces centered at T. At hv = 7.5 eV the shape of the 
inner Fermi surfaces strongly changes and becomes very 
complicated, which is consistent with a bulk state [see 
Fig. S4 of SM] . The size of the energy gap is estimated 
as 230 meV as shown in the 3D map for hv = 10 eV 
[Fig. 1(d)]. Such a large energy gap is beneficial for a 
high stability of the spin current conductance at room 
temperature. Figure 1(d) illustrates a high anisotropy of 
the bulk valence band: the valence band maximum in the 
MfM line is at E B = 490 meV (k\\ = ±0.3 A" 1 ), while 
in the KTK line it is deeper in energy. 

Next we discuss the surface Dirac cone of PbBi2Te4 in 
more detail by comparing it with other 3D TIs. The 
Dirac cone dispersion along PM and TK is shown in 
Fig. 2(a) for PbBi 2 Te 4 and for the well studied TIs Bi 2 Se 3 
and TlBiSe2- Close to Ep the Dirac cone energy dis- 
persion in PbBi2Tc4 is as steep as in the other materi- 
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FIG. 2: (Color online) (a) Experimental energy dispersion 
curves of surface Dirac cones along FM and FK lines in 
PbBi2Fe4 compared to Bi2Se3 and TlBiSe2. Fhe curves are 
obtained from the intensity maxima of momentum distribu- 
tion curves. Fhe energy is relative to the Dirac point, (b) 
Constant energy contours for PbBi2Fe4. (c) Estimated car- 
rier densities S{E)/A-k 2 for three materials as a function of 
energy with respect to Dirac point. S(E) is the area of the 
constant energy contour at energy E. 



als, but it is apparently less steep near the Dirac point. 
In PbBi2Tc4, the group velocity at E-p is estimated as 
3.9 x 10 5 m/s, while near the Dirac point, it is much lower 
(1.4 x 10 5 m/s) than in Bi 2 Se 3 (2.9 x 10 5 m/s) and in 
TlBiSc 2 (3.9 x 10 5 m/s) [jjfa. The sizes of the iso-energy 
contours in the bulk energy gap from the Dirac point to 
200 meV are much larger for PbBi2Tc4 than for the other 
two materials, as shown in Fig .2(b). The estimated topo- 
logical surface carrier density in PbBi2Tc4 obtained from 
the area of constant energy contour S(E)/4ir 2 is much 
larger than in the other two materials [Fig. 2(c)]. 

The theoretical bulk and surface band structures of 
PbBi2Te4 are shown in Figs. 3(a) and 3(b), respec- 
tively. The calculations were performed with the VASP 
code 2j| 3(J with optimized internal lattice parameters. 
In principle, three-dimensional materials with inversion 
symmetry are classified with four Z2 topological invari- 
ants i>o; (^1^2 ^3), which can be determined by the parity 
£m(Ti) of occupied bands at eight time-reversal invari- 
ant momenta (TRIM) r i= ( ni „ 2 „ 3 ) = (nibi + n2b2 + 
n^hs)/2, where bi, b2, are primitive reciprocal lattice 
vectors, and n 3 ; = or 1 [l|, Q. The Z2 invariants are de- 



termined by the equations (— l)' y ° = Yi $i an d (— l)' yfc = 

i=l 
N 

II <5i=(nm 2ra3 ), where 8i = J] 6™( r i) H- For 

n fc = l;n J - 7 £ fc =0,l m=l 

rhombohedral lattice of PbBi 2 Te 4 the TRIMs are T, Z, 
and three equivalent F as well as F points [see SM, 
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FIG. 3: (Color online) (a) Bulk band structure of PbBi 2 Te 4 
calculated along high symmetry directions of the Brillouin 
zone; colors show the weight of the states at Bi (blue), Te 
(red), and Pb (green) atoms, respectively. Signs of Si = ±1 
at the TRIM are also shown, (b) Calculated surface (solid 
line) and bulk continuum states (shaded) of PbBi2Te4 along 
K-F-M line. The inset shows the surface Dirac cones in the 
limited momentum space for f-M and F-K lines. 



Fig. Sl(b)]. The previous study confirmed that this com- 
pound is a strong topological insulator with the principal 
topological invariant uq = 1 [23]. Here we analyze the 
(1/1^2^3) invariants. Interestingly, the parity inversion of 
bulk bands occurs at the Z point for PbBi2Tc4 [Fig. 3(a)], 
which leads to So invariants 1; (111). This is in con- 
trast to the case of binary chalcogenides Bi 2 X 3 (X=Se, 
Te), where the parity inversion takes place at the T point 
with Z 2 invariant 1; (000) 0. Note that PbBi 2 Te 4 is the 
first case among the experimentally established topolog- 
ical insulators with Z2 invariant 1; (111) possessing a 
single Dirac cone surface state. It is, thus, distinguished 
from the Bii-^Sb^ alloy with the same Z2 invariant (3l| 
but with 5 or 3 pairs of surface states crossing the Fermi 
energy @, I n ^h-xSh x , owing to nonzero invariants 
(j/1^2^3), a one-dimensional (ID) topologically protected 
state can exist at the dislocation core [32j . In the case of 
the layered crystal the bulk dislocations can hardly exist, 
but other types of ID TI states, such as edge states in 
thin films or ID states at step edges are possible. 
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Finally, we compare the experimental ARPES results 
with the theoretical band structures, see Fig. 3(b). The 
calculated band structure well reproduces the bulk con- 
duction band minimum located at T point at 200 meV 
above the Dirac point. Although the bulk valence band 
maximum is higher in theory than in experiment, the 
theory well reproduces its location in k-space, that is, 
it appears around l/3fM (ku ~ 0.3 A -1 ). Anisotropic 
features along different symmetry lines can be recognized 
above f 00 meV relative to the Dirac point, as depicted in 
the inset. In addition, the parabolic surface state appears 
at ~ 280 meV above the Dirac point in the conduction 
band gap, which is consistent with the present observa- 
tion [see SM, Fig.S3]. 

Our conclusion led by the present experiment is 
twofold: (i) PbBi2Te4 is proved to be a three-dimensional 
topological insulator with the energy gap of 230 meV and 
a single Dirac cone at the f point, (ii) The size of the 
Fermi surface contour in the bulk energy gap is signifi- 
cantly larger than in the other presently known 3D topo- 
logical insulators, whereby the highest carrier density of 
the known topological surface states is achieved. These 
novel findings pave a way for the efficient control of the 
group velocity with sufficiently large spin current density 
by tuning the chemical potential in the bulk energy gap. 
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